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Abstract 

We report methods for correcting the photoluminescence emission and excitation 
spectra of highly absorbing samples for re-absorption and inner filter effects. We 
derive the general form of the correction, and investigate various methods for deter- 
mining the parameters. Additionally, the correction methods are tested with highly 
absorbing fluorescein and melanin (broadband absorption) solutions; the expected 
linear relationships between absorption and emission are recovered upon applica- 
tion of the correction, indicating that the methods are valid. These procedures allow 
accurate quantitative analysis of the emission of low quantum yield samples (such 
as melanin) at concentrations where absorption is significant. 
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1 Introduction 



Melanins are an important class of dark biological pigment known to be re- 
sponsible for photoprotection of the skin, hair and eyes of many different 
species [1,2,3]. Of the two types found in humans (eumelanin: black, and 
pheomelanin: red-brown) eumelanin is the most common [4]. Paradoxically, 
it is also believed that melanins are directly involved in the UVA-induced 
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photochemical processes that lead to DNA damage in skin cells [5,6,7]. The 
biological structure-property-function relationships of melanins and their con- 
nections to melanoma formation have become the foci of intense scientific 
interest and extensive study. Despite this activity, there remain many unan- 
swered questions regarding its structure and properties [8,9,10]. This fact is 
in part due to melanin's strong, broadband and highly non-linear absorbance 
(Figure 4), very low quantum yield [11] and insolubility [1,12] which make 
spectroscopic analysis challenging. Emission levels are so low that they must 
be measured at concentrations where re-absorption and inner filter effects (at- 
tenuation of the probe beam) strongly distort the resulting spectra (Figures 
3 and 5). It is clear that a procedure to accurately and easily correct for 
these effects is required. This would allow future spectroscopic measurements 
of eumelanin (and other strongly absorbing species) to be analysed accurately 
and quantitatively. Re-absorption correction methods have been previously 
applied to melanin [13]; the work reported here is a far more detailed analysis 
of the assumptions implicit in these methods, and of the effectiveness of the 
procedures. 



2 Calculations 

2. 1 General Form of the Correction 

The following derivation is for a collimated excitation beam incident upon a 
square cross-section cuvette (containing the solution of interest), with collec- 
tion at 90° with respect to excitation (and in the same horizontal plane), as 
shown in Figure 1. This is the most common geometry used, although this 
correction method could easily be apphed to other geometries. We assume a 
small excitation volume in the cuvette (defined by the emission and excitation 
slit widths) so that we can apply the Beer Lambert Law [14] directly. From 
Figure 1, for the incoming and outgoing light, we have: 



where Ii to J4 are the intensity of the beam at the locations indicated in Figure 
1, di and c/2 are the path lengths into and out of the cuvette, and cti = a{Xi) 
and a2 = a{X2) are the absorption coefficients along those paths (dependant 
upon the excitation and detection wavelengths, Ai and A2 respectively). 

Let the subscript m denote intensities affected by re-absorption and inner filter 
effects ('measured' intensities), and the subscript c denote intensities corrected 



l2 = he 
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for these effects. We can then spht the intensities I3 and I4 into components 
as follows: 



/3 = Fe {I2, Ai, A2) + {I2, Ai, A2) + Rc {I2, Ai, A2) (3) 
h = Fm (/i, Ai, A2) + (Ji, Ai, A2) + Rm (/i, Ai, A2) (4) 

where F is the component due to fluorescence of the sample (which we seek to 
recover), B is the background fluorescence (measured from the solvent alone) 
and R is the component due to Raman scattering from the solvent [15]. Ai 
and A2 are the excitation and emission wavelengths respectively. In normal 
situations the fluorescence, Raman scattering and background components 
are all linearly dependant upon the illuminating intensity [14,15], so we can 
assume: 

F,{bI)=bF,{I) (5) 
B,{hI)=hB,{I) (6) 
R,{hI)=hR,{I) (7) 

where h is an arbitrary constant. Combining equations 1 to 7 we find: 

Fe(Ji) = e"i'^^+"^'^^ [F^{h) + BM + i?^(/i)] - [B,{h) + R,{h)] (8) 

In a more convenient form: 

Fcorrected F,irfi£(isured Fijdckground (9) 

where 

fc = /c(Ai,A2) -e"i'='^+"2'^2 (10) 

and Fmeasured (the cmissiou we would actually measure) and Fbackground (the 
background emission we would measure from the solvent alone) are given by 
the bracketed expressions in equation 8. Fcorrected is the quantity we seek - 
the emission due to the fluorophore alone, corrected for inner filter and re- 
absorption effects. Using standard propogation of random, independant errors 
[16], the uncertainty in the corrected photoluminescence will then be given by: 



'^Fcorrected — kFmeasuredy CHiAdi + Q;2^<^2 (H) 

where we have assumed that the uncertainty in di and 0^2 dominates (this was 
found to be the case for our system). 
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cui and 0:2 can be determined from a previously measured absorption spec- 
trum. For a photoluminescence (PL) measurement ai will be a constant (the 
absorption coefficient at the excitation wavelength), whereas a2 will be wave- 
length dependant. For a photoluminescence (PLE) measurement, a2 will be 
constant (the absorption coefficient at the detection wavelength) and ai will 
be wavelength dependant. The only remaining parameters, di and ^2, can be 
determined in a variety of ways, as follows. 

2.2 Method 1: Direct Measurement 

By assuming that the small excitation volume is in the centre of the cuvette, 
di and ^2 can simply be calculated from the known dimensions of the system 
(as shown in Figure 1). For a convential spectroscopic system this assumption 
will be valid, and the geometry well known. 

2.3 Method 2: Raman Peak Attenuation 

If direct measurement of the system is not possible (in thin films for exam- 
ple, or another unique geometry) then we seek another measure of the re- 
absorption. The attenuation of the signal associated with Raman scattering 
of the probe by the solvent molecules [15] (the 'Raman scattering peak') pro- 
vides an alternative method of determining the necessary correction factor. 
From equation 8 we see that 

A;(Ai, \2)Rm{h, Ai, A2) = Rc{h, Ai, A2) (12) 

since only can respond to a change in R^, and i?c = must imply R^ = 0. 
Combining with equation 10 for k we find: 

(0:2) (-Rm) (0:2) ^^"^'^^ 

If excitation at wavelength Xir gives a Raman scattering peak at wavelength 
A2ij we have: 

\Ol[X2R) J \Rm[X2R) J \Oi{X2R) J 

By exciting the solution at a variety of wavelengths we can plot the above 
function such that di and d2 are the gradient and y-intercept respectively. Rc 
and Rm can be determined by fitting the PL spectra to multiple Gaussians, 
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(13) 



as shown in Figure 2. Rc and R,n arc the amphtudes of the Gaussian Raman 
scattering peak in the background and fluorophore spectra respectively. The 
uncertainty in di and d2 can be determined from the hnear regression, and 
used in equation 11. 

This method theoreticaUy allows correction of the emission spectra even if the 
sample geometry is unknown (such as in thin films), as long as the Raman 
peak is clearly visible. Even in the well defined square cuvette geometry, this 
method potentially allows for the fact that the illumination volume may not 
be directly in the centre of the cuvette. 



2.4 Method 3: Raman Peak Attenuation Approximation 



If it were not possible to apply Method 2 (as shall be discussed in the results 
section), then the following approximations may allow correction of unique 
system geometries. If re-absorption dominates over inner-filter effects (0:2^2 ^ 
ai^i) we can make the approximation: 

k = e'^'^^^ii (15) 



where dg// is an effective path length. This may be the case either because 
of the geometry of the system, or because of the absorbance profile of the 
sample. Following the same procedure as for Method 2 (solving for d^ff at the 
excitation wavelength \ir that gives a Raman scattering peak at wavelength 
\2r) we find 

, _ ( 1 \ , ( Rc{\2b) \ 



(16) 



Hence the full correction is: 

Fcorrected{^2) — ©Xp ( — 7T — 4" In ( '- 1 ] F^f,asured{^2) ~ -^6ocfcgrounci('^2) (17) 




with uncertainty: 

«(A2)A:(A2)F„,(A2) 



AF,(A2) 



2R 



( ARrX\2R) 
\ [ Rc{X2r) 



+ 



( ARrn{X2B) 
V -Rm(A2_R) 



(18) 



assuming random, independant errors, dominated by the uncertainty in R^ 
and Rm [16]. Rc and Rm can be determined as for Method 2. 
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Note that although equation 15 no longer explicitly includes the excitation 
wavelength (as the correction factor did in Method 1), d^ff will be dependant 
upon Ai due to the non-linearity of the absorption profile. Hence o?e// must 
be recalculated for PL spectra measured at different excitation wavelengths 
unless the absorption coefficients at these excitation wavelengths (ai) are very 
similar, or 0:2^2 is so much larger than aidi that the difference is negligible. 

This also means that each point in a PLE spectrum will have a different d^ff. 
Since this is impractical to calculate (unless the absorption profile is constant 
at all excitation wavelengths) , we can instead manipulate the geometry of the 
system so that inner filter effects are more significant than re-absorption. This 
can be achieved with a rectangular cuvette, orientated so that the excitation 
path length (di) is significantly longer than the emission path length (^2). We 
can then make the alternate approximation: 



(19) 



givmg: 



d, 



eff 



In 



(20) 



Hence the full correction for PLE spectra is: 



Fcorrected{^l) — GXp { —J\ — T" In f r>'^\ \ ^\ J J Pmeasured{^l) ~ -^6acfcgrow('^l) (21) 



with uncertainty: 



AFc(Ai) 



a(Ai)A:(Ai)F^(Ai) 
"(Alii) 



'A/?e(A 



2R 



(^Rrn{X2R[ 

\ \ Rc{\2r) ) \ Rm{\2R) , 



(22) 



once again assuming random, independant errors, dominated by the uncer- 
tainty in Rf. and [16]. 



Note that since the Raman scattering peak is not clearly visible in PLE spec- 
tra, the effective path length [dejf) must be calculated from the appropriate 
PL spectrum (with the cuvette orientated so that inner filter effects domi- 
nate). Any excitation wavelength is sufficient to determine (ie//; as long as the 
Raman peak is clearly visible. 
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3 Experimental 



3.1 Sample Preparation 

Synthetic eumelanin (dopamelanin) derived from the non-enzymatic oxida- 
tion of tyrosine was purchased from Sigma Aldrich, and used without further 
purification. Eumelanin solutions were prepared at a range of concentrations 
(0.001% to 0.005%) by weight macromolecule in high purity 18.2M(] MilhQ 
de-ionised water. To aid solubility, the pH of the solutions was adjusted us- 
ing 0.01 M NaOH to ~ 11.5, and the solutions gently heated with stirring. 
Under such conditions, pale brown, apparently continuous eumelanin disper- 
sions were produced. Fluorescein was purchased from Sigma Aldrich and used 
without further purification to prepare standard solutions at ten different con- 
centrations varying from 1.2 x 10^'^% to 5 x 10^^% by weight in O.IM NaOH 
solution (18.2Mr2 MilliQ de-ionised water). Fluorescein and eumelanin con- 
centrations were chosen so as to maintain absorbance levels within the range 
of the spectrometer. 



3.2 Absorption Spectrometry 

Absorption spectra between 200nm and 800nm were recorded for the syn- 
thetic eumelanin and fluorescein solutions using a Perkin Elmer Lambda 40 
spectrophotometer. An integration of 2nm, scan speed of 240nm/min and slit 
width of 3nm bandpass were used. Spectra were collected using a quartz 1cm 
square cuvette. Solvent scans (obtained under identical conditions) were used 
for background correction. 



3.3 Photoluminescence Emission Spectrometry 

Photoluminescence emission spectra for the eumelanin and fluorescein solu- 
tions were recorded for all concentrations using a Jobin Yvon FluoroMax 3 
Fluorimctcr. Emission scans were performed between 400nm and 700nm using 
an excitation wavelength of 380nm for the eumelanin samples, and 490nm for 
the fluorescein samples. A band pass of 3nm and an integration of 0.3s were 
used. Background scans were performed under identical instrumental condi- 
tions using the relevant solvents. Spectra were collected using a quartz 1cm 
square cuvette. Spectra were automatically corrected to account for differences 
in pump beam power at different excitation wavelengths using a reference 
beam. 
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3-4 Photoluminescence Excitation Spectrometry 

Photoluminescence excitation spectra for the eumelanin solutions were recorded 
using a Jobin Yvon FluoroMax 3 Fluorimeter. Excitation scans were per- 
formed between 300nm and 465nm using a detection wavelength of 485nm. A 
band pass of 3nm and an integration of 0.3s were used. Background scans were 
performed under identical instrumental conditions using the relevant solvents. 
A rectangular cuvette orientated such that di — 0.5cm, ^2 = 0.2cm was used 
to reduce re- absorption effects in comparison to inner filter effects. Spectra 
were automatically corrected to account for differences in pump beam power 
at different excitation wavelengths using a reference beam. 



4 Results and discussion 

4-1 Method 1: Direct Measurement 

The uncorrected PL spectrum of eumelanin for three different concentrations 
is shown in Figure 3. Inner filter and re-absorption effects are clearly evident; 
the intensity does not increase linearly with concentration, the peak shifts 
to lower energies at higher concentrations, and the Raman scattering peak 
(at approximately 436nm) is increasingly attenuated as the concentration in- 
creases (such that subtraction of the background leaves 'holes' in the spectra). 
The shift in peak position is due to the nonlinearity of eumelanin absorbance 
(Figure 4). 

Figure 5 shows the same PL spectra corrected using Method 1 as outhned 
above, with measured path lengths di — d2 — 0.50cm. The peaks now align 
very closely (to within 2nm), the Raman peak is completely removed by sub- 
traction of the background, and the peak intensity is now linear with concen- 
tration, as shown in Figure 6. 

Figure 7 shows the integrated PL against the absorbance at the excitation 
wavelength for various concentrations of eumelanin. Open circles indicate 
the uncorrected data, and closed squares show the data after the correction 
(Method 1) has been applied. The eumelanin PL is heavily attenuated at 
all but the lowest concentration, but the expected linear relationship is fully 
recovered upon application of the correction. 

As a standard fluorophore, fluorescein was used as a further test of the correc- 
tion method (refer to Figure 8). Due to the very narrow absorbance range of 
fluorescein, unlike eumelanin, inner fllter effects are far more signiflcant than 
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re-absorption. Also unlike eumelanin. the very high quantum yield of fluo- 
rescein means that very low concentrations give measureable emission. Five 
concentrations were used in the low concentration limit [17,18], and five at 
higher concentrations as a test of the method. The emission was attenuated 
as expected at these higher concentrations, but the expected linear relation- 
ship was fully recovered upon application of the correction (refer to Figure 
8). Note that since fluorescein has a very sharp absorption peak (Figure 9), 
exciting at this peak (as wc did here) will make inner fllter effects far more 
signiflcant than re-absorption (aidi ^ a2d2)- 

These measurements indicate that this correction method is valid, both for 
correction of re- absorption, and inner fllter effects (or any combination of the 
two). 

4-2 Method 2: Raman Peak Attenuation 

Figure 10 shows the plot to determine di (as the gradient) and ^2 (as the 
intercept) according to equation 14. The three different solutions (of different 
concentrations) were measured in the same system, and so should all give 
identical values for di and d2 (and hence fall on the same line). Figure 10 
shows that this is clearly not the case; although each data set appears to form 
a straight line, they are all signiflcantly different, and none of the data sets 
gives close to the predicted values for di and ^2 (as indicated by the dotted 
line). Two of the data sets actually have negative gradients, and hence give 
negative values for di. 

The uncertainty in each data point in Figure 10 is very large, indicating that 
the linear trend in each data set is due to a systematic error in the fltting 
routine. Hence we conclude that although this method of correction is theo- 
retically correct, in practice there is too much freedom in the fltting process to 
determine Rc and Rm (and hence di and (^2) accurately. These parameters may 
be better determined by direct measurement (Method 1). It is possible that 
this method may be viable for a sample with a larger Stokes shift; the Raman 
peak would then be separated further from the sample photoluminescence, 
making the fltting process more accurate. 

4.-3 Method 3: Raman Peak Attenuation Approximation 

Figure 11 shows the PL of the 0.0025% eumelanin solution (excited at 380nm) 
corrected for re- absorption and inner filter effects with Methods 1 and 3, as 
well as the uncorrected data (with background subtracted). Note that since the 
Raman scattering peak has been attenuated, subtraction of the background 
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without prior correction leaves the artifact observed at 436nm. The two meth- 
ods agree well close to the Raman peak, but further away the approximation 
used in Method 3 is less accurate (as would be expected). 

In Figure 12 we see the integrated PL against concentration for the two cor- 
rection methods. Both recover the expected linear relationship, and they agree 
at lower concentrations. At high concentrations they deviate more signifi- 
cantly, indicating that if Method 3 must be used, extremely high concen- 
trations should be avoided. Method 3 still offers a vast improvement over the 
uncorrected data, however, as evidenced by Figure 11. 

The effectiveness of Method 3 could not be investigated with the fluorescein 
solutions since the Raman scattering peak was of much lower intensity than 
the dye photoluminescence and hence not visible. Note that this method could 
be applied to any emission or scattering feature that is detectable in both the 
background and sample spectra (not only the Raman peak). 

4-4 Correction of Photoluminescence Excitation Data 

These correction methods can be equally applied to photoluminescence excita- 
tion (PLE) data. For many materials, including melanins, accurate PLE data 
is critical for clear understanding of the emission behaviour. Figure 13 shows 
photoluminescence excitation data for a eumelanin solution. The uncorrected 
data (with background subtracted) shows significant change when correction 
Method 1 is applied. Method 3 does not agree closely, but still provides a sig- 
nificant improvement over the uncorrected data. This is the first publication 
of PLE data for melanin, and will be expanded upon in future publications. 



5 Conclusion 

Two viable methods to correct for re-absorption and inner filter effects in 
emission and excitation spectra have been proposed and validated: 

(1) Direct Measurement (Method 1): This method is the most accurate and 
easy to apply of those proposed. Upon application of this correction all 
of the expected behaviours are recovered. The geometry of the system 
must be well defined. 

(2) Raman Peak Attenuation Approximation (Method 3): This method can 
be used even if the sample geometry is completely unknown (such as 
in thin films) as long as the Raman scattering peak is visible. It is not 
as accurate as the first method, but gives a reasonable approximation. 
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especially close to the Raman peak. 

Method 2 (exact correction via the Raman peak attenuation) failed in this case 
for technical reasons, but may be possible for situations where the Raman peak 
and fiuorophore emission are better separated. 

With these methods it is possible to recover the emission and excitation spec- 
tra of highly absorbing samples (such as melanin), which can be strongly 
distorted by re-absorption and inner filter effects. These methods also make it 
possible to obtain far more accurate estimates of quantities such as the quan- 
tum yield in samples where high absorbance is unavoidable. These methods 
have made it possible to obtain quantitative photoluminescence and photo- 
luminescence excitation data for melanin, and most significantly the first ac- 
curate measurement of the quantum yield of melanin [11]. This is a critical 
step towards understanding the energy dissipation pathways of this extremely 
important biological molecule. 
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Fig. 1. Geometry of the PL and PLE cuvette (view from above). Light at wavelength 
Ai travels a distance di into the cuvette, where it excites photoluminescence (at a 
wavelength A2). This emission then travels a distance d2 through the solution to 
the detector, /i is the initial intensity of excitation light, which is attenuated to 
intensity I2 at the centre of the cuvette (inner filter eff'ect). I3 is the initial intensity 
of the emission from the excitation volume, which is attenuated to intensity I4 when 
it reaches the detector (re-absorption). 
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Fig. 2. PL (380nm excitation) of 0.0025% eumelanin solution with multiple Gaussian 
fitting (dashed: Gaussian components, solid: data and resultant fit, dotted: Raman 
scattering component fit). 
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Fig. 3. Raw PL emission spectra (pumped at 380nm) for three synthetic eume- 
lanin solutions: 0.005% (solid), 0.0025% (dashed) and 0.001% (dotted) by weight 
concentration, and solvent background (dot-dash). 
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Fig. 4. Absorption spectra of synthetic eumelanin solutions at three concentrations: 
0.005% (dotted), 0.0025% (dashed) and 0.001% (sohd) by weight concentration. 
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Fig. 5. PL emission spectra corrected using Method 1 (pumped at 380 nm) for 
three synthetic eumelanin solutions: 0.005% (solid), 0.0025% (dashed) and 0.001% 
(dotted) by weight concentration. Errors were calculated according to equation 11. 
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Fig. 6. PL emission peak intensity (corrected using Method 1) vs. concentration 
for three synthetic eumelanin solutions: 0.005%, 0.0025% and 0.001% by weight 
concentration. The samples were pumped at 380 nm. Intensity errors were calculated 
according to equation 11. 




Fig. 7. Integrated PL emission vs. absorption coefficient at the excitation wavelength 
(380 nm) for five melanin solutions (0.001% to 0.005%): uncorrected data (open 
circles) and corrected with Method 1 (solid squares with linear regression). Errors 
were calculated according to equation 11. 
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Fig. 8. Integrated PL emission vs. absorption coefficient at the excitation wavelength 
(490 nm) for ten fluorescein solutions (1.2 x 10 — 4% to 5 x 10~^% by weight): 
uncorrected data (open circles) and corrected with Method 1 (solid squares with 
linear regression). Errors were calculated according to equation 11. 
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Fig. 9. Absorption spectrum of fluorescein solution (concentration by weight: 
2 X 10-5%). 



17 




Fig. 10. Plot to determine di and d2, as described in the text. Three solutions 
of different concentrations were measured: 0.001% (open circles), 0.0025% (filled 
squares) and 0.005% (open triangles). The dotted line indicates the predicted trend, 
calculated from the known system geometry. Errors were calculated according to 
equation 11 as described in the text. 
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Fig. 11. PL of 0.0025% eumelanin solution (380nm excitation). Solid: corrected 
using Method 1, dashed: corrected using Method 3, dotted: no correction applied 
(background subtracted giving the artifact at 436 nm). 
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Fig. 12. Integrated PL of eumelanin solutions (380nm excitation) corrected with 
Method 1 (closed squares), Method 3 (open triangles) and uncorrected with back- 
ground subtracted (open circles). Linear regressions for each data set (shown) were 
constrained to pass through the origin. Although at high concentrations Method 3 
does not agree with Method 1 to within the uncertainty, it still provides a signif- 
icant improvement over the uncorrected data. Errors were calculated according to 
equations 11 and 18. 
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Fig. 13. Photoluminescence Excitation of eumelanin solution (0.001% by weight). 
Uncorrected, with background subtracted (dotted) , corrected with Method 1 (sohd) 
and corrected with Method 3 (dashed) . The Raman peak height ratio for Method 3 
was calculated from a PL spectrum excited at 380nm in the same geometry (rect- 
angular cuvette with inner filter effects dominating). 
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